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E APPLICATION OF GROUP THEORY 
TO THE NORMAL VIBRATIONS: 
OF A CUBIC CRYSTAL* 


ABSTRACT 


HE use that can be made of group theoretical considera- 
tions in analyzing the normal vibrations of a crystal is 
lustrated in the case of a simple cubic lattice. The group 
nder which the system is invariant can be made finite by 
he use of periodic boundary conditions. The invariant sub- 
roup containing the translations only can be completely 
teduced on the basis of plane waves of arbitrary polarization. 
hese can be divided into sets of 144 members, and each set 
ives rise to a representation of the whole group of the crys- 
al. In a number of special cases these representations can 
be reduced with respect to various sub-groups of rotations 
to such an extent that it is possible to determine the direc- 
tions of the normal vibrations. In such cases the frequencies 
can be directly expressed in terms of the force constants. 

The group theoretical treatment makes it easy to distin- 
guish those properties of the elastic spectrum which depend 
only on the symmetry of the lattice, from those which de- 
pend on special properties of the force constants. 

* * * 

The principal applications of group theory to mathemati- 
cal physics are based on the invariance of the differential 
equations of physics to groups of transformations. The vari- 
ous solutions of the equations can be used as the elements 


*By William Vermillion Houston, Ph.D., Professor of Physics and President of 
the Rice Institute. 
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on which a representation of the group can be based, and 
from the properties of the group and the theory of group 
representations, it is possible to draw some conclusions con- 
cerning the solutions themselves. I propose to illustrate here 
this use of group theory by considering in some detail its 
application to the problem of the normal vibrations of a 
simple cubic crystal lattice. 

The problem of normal vibrations is of physical impor- 
tance in a number of ways. It provides the basis for com- 
putation of the specific heat of the crystal, and for the effect 
of temperature on x-ray diffraction. It is also of importance 
in the theory of the electrical conductivity of metals and in 
the theory of residual rays. The first approach to such prob- 
lems was made by Debye in 1913, when he developed his 
theory of specific heats. In this work he ignored all of the 
characteristically crystalline properties of the solid and 
treated it as an elastic continuum. About the same time 
Born and von Karman showed how the matter could be 
treated in terms of normal vibrations, and recently the ap- 
plication, of group theory has made it possible to draw some 
general conclusions concerning the elastic spectrum of a crys- 
tal with a minimum of calculation. Although these conclu- 
sions give one a good view of the general situation, they usu- 
ally fall just short of giving all of the information which is 
desired and need to be supplemented by some numerical 
calculation. Several years ago P. C. Fine, at the California 
Institute of Technology, carried through a numerical evalua- 
tion of the frequencies of the normal vibrations of a body- 
centered cubic crystal. As far as I know this is the first 
time that such a calculation has been carried to an explicit 
conclusion, and it provides an interesting example of a fre- 
quency distribution. However, the methods which I shall 
describe here are based on the group theoretical work of 


7 ert) 7 ae 
iy ‘pt eS vibe : p i” 74D Per 4) beeiantie 


th qviari s wor 7 » naamcA Boe. fae iol mf 


* . P ~*~ 
att ae bie f) ', f. “ St, d mary ie ae VOR at Lee tri > ~ 2 i. 
. ; ‘ | | ; +e Se a) 
eH OF SLe8O9 34 then Bo yeas UTR MOTI Q 
oS Ve etuvioe jreale ocly gANtYootOn pi ile tt mi 
: ; 


“..*® 
64 Loni * ania 


s radio cnrem's, word? bere: 28 
© Ya trpot cits lenein-edt Jo yeti: 
dy" en aly wegen Eine ta) A ‘tawge 
at 0 OA riko Aga? ase Tonwiialag, & 
0 sfoyrumam jenireone ort ei babligy 


sare fe aed saan tale ove 
ott 0 bon 
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Wigner,’ Seitz,” and others, and lead to only those results 
that can be obtained from the symmetry of the problem 
and are independent of the nature of the forces. 


THE NORMAL VIBRATIONS 


For purposes of illustration, and to avord complications, 
consider a simple cubic lattice. In its equilibrium configura- 
tion, this consists of a particle of mass m at every point 
whose coordinates, in a Cartesian system, are integers. The 
masses are connected together in such a way that the force 
- acting on one of them is a linear function of the displace- 
ments of all of the points, itself included. If r,,,,, is the dis- 
placement from its equilibrium position of the particle whose 
rest position has the coordinates (/,m,n), then the differen- 
tial equations of motion are 


d* m,n m,n 
moe —— = > Fete ans (1) 
A,u, 0 


The constants fates depend only on (J—X), (m—y), (n—»), 
and get smaller as (J—\)?+(m—y)?+(n—v)? increases. In 
order to avoid complications due to the boundaries, it is 
customary to apply cyclic boundary conditions. By these 
it is required that 


Pimn = 1 14+aG,m46,n4cG (2) 


where G is a large integer that determines the size of the 
cycle, and a, b, and c are arbitrary integers. This type of 
boundary condition eliminates immediately all phenomena 
connected with the crystal surfaces. A model such as this 
will describe only pure volume phenomena, but it is believed 


1E, Wigner, “Gruppentheorie,” Vieweg u. Sohn. 

2F. Seitz, “A Matrix-algebraic Development of the Crystallographic Groups,” 
Zs. f. Kristallographie, vol. 88, p. 433, 1934; vol. 90, p. 289, 1935; vol. 91, p. 336, 
1936; vol. 92, p. 100, 1937. F. Seitz, “On the Reduction of Space Groups,” Annals 
of Mathematics, vol. 37, p. 17, 1936. : 
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to include all of the elements essential to a discussion of the 
normal vibrations. Since it is linear, this model does not 
provide for the treatment of thermal conductivity or thermal 
expansion. These phenomena are associated with departures 
from linearity in the differential equations of motion. The 
boundary conditions reduce the number of effectively dif- 
ferent lattice points to G*, and G must be taken so large 
that the results are essentially independent of its size. 

The existence of normal coordinates for a set of equations 
of this form is given by the general theory. These coordinates 
are linear combinations of the displacements, 


Ps — >> Cie Fimn (3) 
such that 
dp, oe 
a= — 8 Ps (4) 


The solution of these equations can be written down at once, 
and then by inversion of equations (3) 


Timn a oe: Pie P; (S) 


This gives the final solution of the problem. The principal 
point is then actually to find the normal coordinates, and 
the various coefficients C;,,,. By the application of group 
theoretical considerations this can be reduced to the solu- 
tion of a cubic equation, when the lattice is simple cubic, 
and to an equation of order equal to three times the number 
of masses in a unit cell for other kinds. In special cases 
these equations can be further reduced by additional con- 
siderations based on the rotational symmetries of the crystal. 
This is the object of the group theoretical discussion given 
here. 

Vhe differentia! equations (1), and consequently (4), are 
invariant to any translation or rotation of the crystal lattice 
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at transforms the various equilibrium positions into them- 
ves. The content of this statement is simply that the co- 
cients 2 Babe depend only upon the differences of the in- 
ces. The translations are 


l+Il-+-a, m-+>-m-+b, n>n+e. T: 


he rotations are, for the simple cubic lattice, the members 

the rotation group O to which a simple cube is invariant. 
he differential equation is also invariant under the inver- 
on, or reflection through the origin, in which +—/, m+—™m, 
-+—n. The normal coordinates p, constitute a basis for 
he establishment of a representation of this whole group 
ecause the differential equation for a normal vibration is 
nvariant under all these transformations. Furthermore, 
ince each differential equation for a normal vibration in- 
ludes only one frequency, the effect of any of the transfor- 
nations of this group will be to transform one normal vibra- 
ion to another, or linear combination of others, which 
ibrates with the same frequency. Hence those normal vi- 
rations whose frequencies are identical constitute a basis 
or a representation of the whole group of the crystal. For 
his reason it is possible to draw some conclusions about the 
ormal vibrations from the properties of the group and its 
rreducible representations. 


THE GROUP OF A CUBIC CRYSTAL 


In discussing the properties of this group I shall go into 
onsiderable detail in giving illustrations which have been 
istructive to me. To those who are well versed in group 
heory these things may seem superfluous, but to me a proof 
; more convincing when accompanied by an illustration. 

The group of translations constitutes an invariant sub- 
roup of the whole group. This means that any rotation, 
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followed by a translation, and then by the inverse of the 
original rotation, is just equivalent to a simple translation. 
This invariant sub-group of translations is an Abelian group. 
All its members commute with one another. It makes no 
difference in what order they are carried out. 

The group of rotations of the crystal is just the group of 
rotations about its center under which a simple cube is in- 
variant. These are represented in Figure 1 as the trans- 
formation of the coordinates x, y, z. The association of these 
matrices with the rotations of the cube can be seen from the 
figure. The first element of the group is the identity opera- 
tion. This leaves everything unchanged. The second ele- 
ment, Ds, represents a rotation of 180° about the x-axis. . 
It changes the sign of the y and z coordinates and leaves the 
x-coordinate unchanged. It is clear from this that a second 
application of this operation will restore the cube to its 
original configuration; so the element is of order 2. This is 
indicated by the subscript. The superscript merely indicates 
that this is the first of the elements of order 2, and the 
letter D indicates that it is a pure rotation. The element 
R; can be regarded as a rotation about the x-axis of 180° 
followed by the change in sign of all three coordinates, or 
the reflection through the origin. This is the justification 
for the notation. The element R:=/D;. A simpler view of ' 
this element is that it is a reflection in the y-z-plane. 

As another example take the element D;. This is a rota- 
tion through an angle 27/3 about the axis symmetrically 
between x, y, and z. By this operation x’ =z, y’=x, z’=y. 
This is obviously an operation of order 3. (There are also 
operations of order 4 which consist of rotations through 90° 
about axes parallel to the coordinate axes.) . 

The 24 elements which represent pure rotations constitute 
an invariant sub-group of the whole. When this invariant 
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sub-group is augmented by the element /, the inversion 
through the origin, there results the complete group O* of 
48 elements. 

The properties of a group are all given by its multiplica- 
tion table. Figure 2 gives the multiplication table for the 
group O. This is, of course, just 1/4 of the table for the 
whole group but the rest can be easily inferred from that 
much. This table is the real definition of the group. The 
association of the elements with the rotation of a cube, or 
with the transformation of coordinates, is merely a repre- 
sentation or an illustration of the group. Nevertheless, it is 
the representations that are of interest in applications to 
physics. 

Figure 3 shows some of the properties of this group based 
entirely on its multiplication table. There is first the divi- 
sion into classes. If an element 4 is multiplied on the right 
by another element B and on the left by the inverse of B, 
the result.is another element in the same class. Such a _ 
process is called the transformation of the element 4 by the 
element B. The transformation of any element 4 of a group 
by any other element of the group leads to an element 
which belongs to the same class as 4. Thus the transforma- 
tion of the element Di by D; leads to Dj, which is a member 
of the class C, as was Dj. 

Any group of elements which constitutes a group by itself 
is called a sub-group. The unit’ element is a sub-group, but 
a trivial one. The unit element and the element D; consti- 
tute a sub-group whose multiplication table is given in the 
figure. An invariant sub-group is a sub-group which con- 
tains only whole classes. Every member of an invariant sub- 
group is transformed into another member of the sub-group 
by any element of the whole group. The group composed 

of the unit element and D; is a sub-group, but it is not an 
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ariant sub-group since it contains only one of the ele- 
nts of the class Ci and that element, D3, is transformed 
0 D; by D;. The group G is the smallest invariant sub- 
up larger than the identity element itself. The next larger 
ariant sub-group is 7, which is composed of G and the 
ments of the class C3. It can also be constructed by multi- 
ing each element of the sub-group G by D; and then 
)?. The group 7 contains three times as many elements 
the group G. 

t is also possible to regard the operation of multiplying 
h element of the sub-group G by D} as a group operation. 
it is called 4 and its square B, the multiplication table 
hown in the figure. Such a group is called a factor group, 


1 is designated by a At the bottom of the figure is shown 


composition series of O*. O is the largest invariant sub- 
up in O*. The whole group O* can be obtained from O 


h 
adding OX /. The factor group o has just two elements; 
it is of the order 2. T is the largest invariant sub-group 


O and Gis also of the order 2. Similarly G is the largest 


rariant sub-group in 7 and, as shown, is of the order 3. 
is the largest invariant sub-group of G. One should note 
1t F is an invariant sub-group of G, but is not an invariant 
»-group of the whole group O*. Transformation of F by 
ments of G leaves it invariant, but transformation of F 
elements outside of G but still in the group O* does not 


ve it invariant. The factor group : is of order 2 and 


ally the factor group % is of order 2. Such a series of 


rariant sub-groups is called a composition series. When 
of the factor groups are of prime order, as in this case, 
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the whole group, O*, is said to be resolvable. Seitz has shown 
that all groups represented by crystals are resolvable. This 
property of resolvability is of importance in the reduction 
of the representation. 

The set of linear transformations, in terms of which this 
group of rotations was specified, constitutes a representation 
of the group. The transformation matrices obey the same 
laws of multiplication as the group elements. It is possible, 
however, to construct many other representations. Suppose, 
for instance, that instead of the coordinates x, Y, z one takes 
all their quadratic combinations x?, y?, z?, xy, xz, yz, and 
applies to them the transformation indicated in the first 
representation. The matrices of these transformations will 
contains six rows and six columns, and will constitute a 
representation of the group. Some examples of these ma- 
trices are given in Figure 4. This shows one representative 
from each class and in the bottom row are three from the 
class JC3. 

A common property of all these matrices is that they con- 
sist of two 3X3 matrices along the diagonal. The upper 
right and lower left squares consist entirely of zeros. This 
fact is described by saying that the representation is partly 
reduced. If it were possible, by using linear combinations 
of these quadratic terms instead of the simple products 


- themselves, to reduce all the matrices to smaller square 


matrices along the diagonal, the representation would be 
further reducible. As a matter of fact, this representation is 
reducible. The matrices in the lower right-hand corner can- 
not be reduced, but those in the upper left corner can be. 
This is shown in Figure 5. The new variables are 


R=x2+y?+2? S=x?+ey?+e” T=x?+ey’+ex? 


and each of the matrices consists of one single term and a 
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Normal Vibrations of a Cubic Crystal 


Ficure 4—A REPRESENTATION OF GROUP O* 
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A REDUCED REPRESENTATION OF O% 
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2 matrix on the diagonal. In this way the original rep- 
entation containing six rows and six columns can be con- 
ered as composed of a one-dimensional, a two-dimen- 
nal, and a three-dimensional irreducible representation. 
is is the limit of possible reduction in this case. 

ach group has a set of irreducible representations. For a 
ite group, the number of such representations is finite. 
yur of the irreducible representations of this group O* have 
en already shown. The original matrices constitute one, 
d the other three are included in the six-row, six-column 
atrices just discussed. 
An important theorem in the theory of representations 
that every representation is composed, in a unique way, 
irreducible representations. For instance the six-row, 
column representation just discussed has been shown to 


- composed of one i-dimensional, one 2-dimensional, and: 


ie 3-dimensional irreducible representation. It might be 
rought that this reduction could be carried out in many 
ferent ways, but such is not the case. A given representa- 
yn can be reduced to only one set of irreducible represen- 
tions. It is largely because of this theorem that the theory 
group representations is useful in treating the equations 
physics. 
Any representation of a group can be described by means 
its group characters. These are the diagonal sums of the 
atrices of the representation. This diagonal sum is the 
me for all the members of a class, so that it is necessary 
give the character merely for each class, and not for each 
ment of the group. In the chart showing the six-dimen- 
nal representation, the identity matrix has the character 
The identity will always have a character equal to the 
mber of dimensions of the representation. The class Ci 


s the character 2 and the other characters are indicated’ 
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on the chart. In particular, the three elements of the class 
JC; are all shown to.have the same character 0. 

When the characters of the representations are given, the 
representation is essentially defined. Matrix transformations 
of the form 4-'B4 lead to what is known as an equivalent 
representation. Such representations are regarded as not 
being essentially distinct from the original even though the 
form of the matrices will be different. Such transformations 
do not change the characters of the representation. For this 
reason, new representations which have the same set of char- 
acters are either identical or equivalent representations. 

Each irreducible representation of a group has a set of 
group characters which uniquely define it. In the chart 
showing the identity representation and the two-dimensional 
representation, it can be seen that all the characters of the 
identity are 1. In the two-dimensional representation the 
character of the unit element is naturally 2; for the class 
C, the character is also 2, as shown in the two examples. 
For classes Cy, Cy, JCs, and JC, the character is zero; for 
class C3 and class JC; it is —1, and for J and JC; the char- 
acter is 2. The specification of this set of numbers, together 
with the classes to which they belong, completely specifies 
this particular irreducible representation. The characters of 
the irreducible representations of O* are shown in Figure 6. 

If a representation is composed of a number of irreducible 
representations, the character of each class is the sum of the 
characters of the various irreducible representations of the 
class. Furthermore, the properties of the characters are such 
that this sum can be made up in one and only one way, cor- 
responding to the unique reduction of the representation. 
In this way it is possible to determine the irreducible repre- 
sentations of which any representation is composed, pro- 
vided the characters are known. 
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GROUP REPRESENTATIONS BASED ON THE 
NORMAL VIBRATIONS 


Let us now return to the vibrations of a cubic crystal and 
attempt to form a representation of the whole group on the 
basis of them. In the first place,,one could use the displace- 
ments of the individual lattice points as the basis of such a 
representation. Each point has three possible components 
of displacement and there are G* points; so the representa- 
tion would be 3G* dimensional. To start to reduce this rep- 
resentation, consider first the invariant sub-group consisting 
of the translations alone. Since the group is Abelian it has 
only one-dimensional irreducible representations ; and since 
every element is of the order G, each matrix must be diagonal 
with Gth roots of unity on the diagonal. The linear combina- 
tions which are the basis for such a reduced representation _ 
are 


2% (kil-Heem Hen) 


P3,kakeayka — Lu Tj,,m.n 
l,mn 


A translation of the lattice obviously multiplies each of these 
quantities by a root of unity, and a G-fold iteration of any 
translation restores the variables to their original values. 
This constitutes a completely reduced representation of the 
invariant sub-group of translations. Each p is characterized 
by an index 7 which can take on the values 1, 2, 3 represent- 
ing the three possible directions of displacement along x, y, z, 
and by the three integers 4, ko, ks. These integers can take 


on values from -¢ to +¢ and the triplets of integers may 


‘be represented by vectors in the ’;, ke, ks space. The ends 
of the possible vectors fill uniformly a cube of edge G. 
If a rotation element of the group O* is applied to one of 
these ;,i.2, it will, in general, transform it to a different 
ky’, ke’, Rs’ and a linear combination of the indices 7. Since 
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he differential equation of motion is invariant to such a 
otation, the time dependence of the transformed displace- 
ient will be the same as that of the original. If such linear 
ombinations of the indices 7 are formed as to give a normal 
oordinate, the rotation will transform it into another nor- 
1al coordinate, since the relative displacements are un- 
hanged by the rotation. The transformed normal coordi- 
ate will also have the same frequency as the original. For 


general vector (k,k2k3) there are 48 directions in which it 
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can be rotated by the members of the group, and so there 
are 48 normal vibrations with the same frequency of vibra- 
tion, and in general there are three frequencies for each 
value of kykok3. This, however, gives no direct information 
_about the frequencies or the form of the vibrations, except 
to show that only 1/48 of the whole cube need be considered 
in computing them. In special cases, however, some such 
information can be obtained, and by making use of the fact 
that the frequencies and the forms of vibration vary con- 
tinuously with &, it is possible to infer some additional in- 
formation. 

Consider first the case where k=0. This is designated by 
I in Figure 7. A transformation of pj) by any of the ele- 
ments of the group O* leaves the value of & invariant and 
affects only the index 7. The group of the vectors ending 
at I’ is the whole cubic group. The representation based on 
these coordinates will be a three-dimensional representation 
and will be just the representation shown for the description 
of the operators. It is based on x, y, and z and is an irredu- 
cible representation. This means that for k=0, vibrations 
in any direction are normal vibrations and they all have the 
same frequency. At this point the three possible frequencies 
of vibration coincide. 

Consider next those vectors k for which kk =k; =0. These 
terminate along the line A in the figure. There is a sub- 
group of the group of rotation to which these wave vectors 
are invariant. The elements and properties of this group 
are shown in Figure 8. The characters of the representation 
of this sub-group that is based on x, y, % are given in the 
first row. Then follow the characters of the irreducible repre- 
sentations of A and finally the composition of the representa- 
tion at hand in terms of the irreducible representation. A, is 

a one-dimensional representation and is based on p;. This 


4, 
rAY 


. “ Mid 
~ 
» ? 
a. a ' 
- ' 
\ 
) 
, 
~ Ses 
a i 7 mS 4 iy : 
F ay 
: Brit 1 ethom a 7 
7 F a 
i te > o> ae 


| wT .qung ofdwe led sda ePieaa 
{ ; L ' q snllnieny § ota m 
: ots b sa 1 seis ey od thw: bas . 
. ? ao band ai 2) page ty 

ry id (he @ icon 26h visa Smt 
¥ Tt ‘i liv lemon perirenad 


- +2 
. 0 wourls Sala onde shy a] 
a _t 

‘ Po 


tthe Ae gh slider ao At 20-eebgs alien: 
din i] saan as ni & ool mw 
4107038V aYWY ? hourw O34 ier age s 


ride sa cnteareysy addy To 8% ae sits eA 5 
adr ni Oevig 40 @ ee beasd a ae: 
~suprt ae iuberd aed7 aeierie vette a@ 2 
au terse Vy) 7 lo mfetieo Leo ait cu hal 
oh ote non 193 sigrs bani eld i 


eat «A, VQ.bre “ail cal Thue satis} i et 
Pi 


7 


3 a 
io - 
* i : 
ad = 
Gel 7 ‘ 
: as er 


Ce een 


ae 


Normal Vibrations of a Cubic Crystal 145 
CHARACTERS OF THE SUB-GROUPS A AND S 
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means that a vibration along the x-axis is a normal vibration 
with a frequency characteristic of it. A; is a two-dimensional . 
representation based on p, and p,. Any linear combination 
of p, and p, is a normal vibration, and the two independent 
linear combinations have the same frequency. A similar 
consideration of the group to which the vectors ending at X 
is invariant shows that at X there are still only two inde- 
pendent frequencies. All along the line A the vibrations have 
the nature of transverse and longitudinal waves as in an elas- 
tic continuum, and there are two independent frequencies. 
Again the vectors that end along the line S are invariant. © 
- to a sub-group S whose properties are shown. This sub- 
group can be completely reduced as shown, and the reduced 
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TWO HOLOMORPHIC SUB-GROUPS OF ©’ 
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presentations are based on p,+p,, pz) py—pz. There are 
us three different frequencies. In this case there are no 
ngitudinal or transverse vibrations, but the normal vibra- 
ns are in the face of the cube and in the x-direction. 
The vector ending at the point R is invariant to all mem- 
ts of the rotation group O*, and so just as at I’ there is 
t one frequency. All of the three possible vibrations at 
is point have the same frequency. Lines A and 7 have 
lomorphic groups (Figure 9). 
Whenever it is possible to find a wave vector k whose group 
ntains more than the identity, it is possible to solve the 
ubic equation which determines the frequencies of vibra- 
ion. Thus over the surfaces of the portion of the & cube 
hat it is necessary to consider, the application of the theory 
f group representations indicates the method of treating 
he cubic equation and provides a basis for classification of 
he roots. In case the representation can be completely re- 
uced, the group theory itself gives the solution of the cubic 
quation. In cases in which there is still a two-dimensional 
epresentation left, the theory of group representation gives 
ne root and the remaining quadratic equation can be easily 
andled by ordinary methods. For those points inside of 
his region, some estimate of the frequencies can be obtained 
y extrapolation from the surface. 
This example which I have developed in some detail is, 
1 itself, not very important. There exists no crystal with 
imple cubic symmetry and the results for such a crystal 
ould be obtained by more direct, if less elegant, means. 
lowever, I think it serves best to illustrate a method which 
an be applied to crystals of all types and permits an imme- 
iate recognition of those properties of the elastic spectrum 
hich are to be attributed to the symmetry of the crystal 
lone, and are independent of the nature of the forces. 


W. V. Houston. 
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